Hydrogen storage in MOF-5: A van der Waals density functional theory study 
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Physisorption of hydrogen molecules in metal-organic frameworks (MOFs) provides a promising way for 
hydrogen storage, in which the van der Waals (vdW) interaction plays an important role but cannot described 
by the density functional theory (DFT). By using the vdW density functional (vdW-DF) method, we investigate 
systematically the binding energies of hydrogen molecules in MOF-5 crystal. We first examine the accuracy of 
this methodology by comparing its results with those from the correlated quantum chemistry methods for several 
fragment models cut out from the crystal. Good comparable accuracy is found. By performing calculations for 
the true crystal structure adsorbing one or multiple H2 in the primitive cell, we show that these fragment models 
which have been focused previously cannot represent well the property of the crystal which cannot, however, 
be dealt with by the quantum chemistry methods. It is found that the binding energy with the organic linker is 
much smaller than with the metal oxide corner, which limits the H2 loading. We show that this can be improved 
significantly (from 5.50 to 10.39 kJ/mol) by replacing the H atoms of the organic linker with F atoms which 
cause extra electrostatic interaction. 



I. INTRODUCTION 

Hydrogen as a substitution for fossil fuel in applications of 
energy carrier has stimulated extensive research because of its 
abundance on earth and especially its clean outcome (water). 
However, a safe and efficient hydrogen storage technique still 
remains as an obstacle in practice, in spite of various chem- 
ical and physical schemes proposed |l-4], in which hydro- 
gen is stored in atomic and molecular forms, respectively. In 
this regards, metal-organic frameworks (MOFs) was found to 
be a promising materials family for hydrogen storage based 
on physical adsorption. So far, a series of MOFs has been 
synthesized and their structural properties and hydrogen load- 
ing under different temperatures and pressures have been in- 
vestigated by powder x-ray diffraction (XRD) and inelastic 
neutron scattering (INS) techniques |l5l Compared with 
other schemes, MOFs have several outstanding advantages, 
such as very large surface/volume ratio, tailorable chemical 
composition and structure, and the regular high-symmetry pe- 
riodic crystal structure. MOF-5 as the first one in the series 
has become a particular focus of MOFs study 1I5I-IT3I1 . It has 
been found experimentally that MOF-5 with the cubic three- 
dimensional extended porous structure can achieve a gravi- 
metric H2 capacity up to 4.5%wt at 78K 0]. However, 
a common problem exists with MOFs for hydrogen storage: 
Around room temperature the H2 binding energy is still not 
large enough. How to increase the adsorption enthalpy to 
above 15 kJ/mol fl4'-fr6] under room temperature has become 
a challenge in practice, which may be possible by tailoring the 
chemical species and atomic structure. To achieve this goal 
we need to first gain a full understanding of H2 interaction 
with MOF-5 on the different adsorption sites in the internal 
space of MOF-5. 

To accurately calculate these interactions is, however, a big 
challenge to theoretical studies because the van der Waals 
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(vdW) interaction plays an important role in the physisorp- 
tion, which is a quantum dispersion force and is related to the 
nonlocal correlation effect. This nonlocal correlation is ab- 
sent in the density functional theory (DFT) and can only be 
described by the many-body perturbation theory or the cor- 
related wavefunction methods. The lowest level of theory 
which can deal with the vdW interaction is the M0ller-Plesset 
second-order perturbation (MP2) theory [17]. More sophisti- 
cated methods include the coupled-cluster method, such as the 
coupled-cluster singles and doubles and noniterative triples 
(CCSD(T)) 1 18]. These correlated quantum chemistry meth- 
ods also face challenges for the study of hydrogen storage in 
MOFs. First, these methods cannot deal with periodic crys- 
tal structures. Consequently, one has to use some fragment 
models cut out from the crystal structure to mimic the differ- 
ent parts of MOFs [T^ . Considering the long-range 
nature of the electrostatic and vdW forces present in the ph- 
ysisorption, whether and to what extent these fragment mod- 
els can represent the behavior of the crystal is a question to 
answer. These correlated methods are also computationally 
very demanding and cannot deal with large fragment models. 
Furthermore, as usually implemented with the Gaussian basis 
functions, these methods have a strong basis set dependence 
and, therefore, the result becomes reliable only when a very 
large basis set is used or some elaborate correction techniques 
are applied 1I21I1 . 

To simply add empirically a 1 /r^ term to conventional DFT 
calculations, as did in the DFT+Dispersion (DFT+D) method 
1 22, 23], may be a solution for studying the H2 adsorption in 
the full crystal structure, but its accuracy and reliability of re- 
sult still needs to be examined On the other hand, 
in the vdW-DF method developed by Dion et at. iH HI the 
vdW interaction can be included perturbatively in a seamless 
fashion into an first-principles Kohn-Sham DFT calculation. 
This is achieved by expanding the correlation energy to the 
second order in terms of a carefully chosen quantity for the 
long-range part of the correlation functional. This method 
has been applied to different separational systems ||2ll|24l] and 
shown reasonable accuracy. For the extended crystal structure 
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of MOFs interacting with H2, whether and to what extent the 
binding energy can be determined by the vdW-DF method has 
not been studied in the Hterature, to our best knowledge. 

In this work, we adopt the vdW-DF method to study sys- 
tematically the binding energy of H2 in MOF-5 crystal for 
two major adsorption sites. We first test the accuracy of 
the methodology by considering the fragment models already 
studied in the literature. The comparison with the MP2, 
CCSD(T), and DFT+D calculations shows a good accuracy 
which is better than the empirical DFT+D method. Calcula- 
tions for the full crystal structure, however, show that these 
fragment models cannot represent well the hydrogen- storage 
property of the crystal due to the long-range interaction of H2 
with the environment. Besides the single H2 binding behavior, 
multiple H2 adsorption is also considered to investigate the ef- 
fect of H2 density. Furthermore, to illustrate the possibility of 
enhancing the average H2 binding energy in MOF-5 crystal, 
we consider the substitution of F atoms for the H atoms of 
the organic linker part. It is found that in this way the H2 
binding energy can be almost doubled due to the extra elec- 
trostatic interaction, indicating a possible way to improve the 
performance of MOFs for hydrogen storage. 



II. THEORY, MODELS, AND COMPUTATIONAL 
DETAILS 

Electrostatic force and vdW force are the sources of the 
physisorption for the hydrogen storage in MOFs. Therefore, 
both of them are important for getting computationally accu- 
rate binding energy of H2. Since the vdW interaction is ab- 
sent in DFT we adopt the recently developed vdW-DF method 
which includes the nonlocal electronic correlation effect into 
the Kohn-Sham DFT calculation in a seamless fashion 1I24I] . 
This is achieved by expanding the correlation energy to the 
second order in terms of a carefully chosen quantity for the 
long-range part of the correlation functional. In this way, 
long-range dispersion effects are included as a perturbation to 
the correlation term of the local or semilocal density approx- 
imation 1I25I1 . This method has been successful in describing 
various groups of materials, ranging from molecules to sur- 
faces and bulk materials [21, 24, 25]. Considering the limita- 
tion of the sophisticated quantum chemistry methods for pe- 
riodic systems, this method may provide a useful tool for the 
study of hydrogen storage with MOF crystals, which has not 
been reported in the literature to our best knowledge. 

The crystal structure of MOF-5 is composed of in- 
organic [0Zn4]^+ octahedral groups linked by 1,4- 
benzenedicarboxylate (BDC) organic linkers, forming a 
highly microporous cubic framework, as shown in Fig. [Ha). 
It has a symmetry of face-centered cubic Fm3m with an 
experimental lattice constant of 25.89 A [5] for the con- 
ventional cell. The conventional cell contains eight formula 
units Zn40(BDC)3 (424 atoms in total) while the primitive 
cell contains only two (106 atoms in total), as shown in 
Fig. [U The large empty hole formed can be an ideal space 
for hydrogen storage where the physisorption takes place. 
Specifically, both the 0Zn4 cluster in the corner and the 
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FIG. 1: (a) The conventional cell (424 atoms) and (b) primitive 
cell (106 atoms) of the MOF-5 crystall which has the symmetry of 
Fm3m and contains eight and two formula units (OZn40(BDC)3), 
respectively. Red, blue, gray, and white denotes O, Zn, C, and H 
atoms, respectively. 



organic linker in the pillar can interact with H2 molecules 
with different interaction strengthes. Here, we consider two 
adsorption sites, one is above the C3 plane of the corner (the 
cup site, see Fig. [3a)) and the other is above the benzene 
ring of the organic linker (linker site, see Fig. [3b)). Previous 
theoretical and experimental studies showed that the cup site 
is the most favored one and will be taken up by hydrogen 
molecules first, then the linker site will be occupied. |6, 31] 

In order to examine the reasonability of fragment models in 
representing the different parts of the full crystal, three frag- 
ment models, corner model, BDC model, and edge model, are 
considered, which are used to represent the metal oxide cor- 
ner, organic linker, and one edge of MOF-5, respectively (see 
Fig. [2]). These models are also used to examine the accuracy 
of the vdW-DF method by comparing its results to the exist- 
ing results from the correlated quantum chemistry methods as 




FIG. 2: Optimized structures of the three fragment models: (a) the 
corner model with a H2 adsorbed on the cup site, i.e., above the C3 
plane indicated by the dashed-line triangle, (b) the BDC model with 
a H2 adsorbed above the benzene ring, and (c) the edge model which 
is made up of the two corner models linked by one BDC linker. 



well as DFT+D. 

Our calculations are carried out by using the pseudopoten- 
tial plane-wave formulism of the vdW-DF theory, as imple- 
mented in the Quantum ESPRESSO packageL26]. The O, Zn, 
C, and H atoms are described by the ultrasoft pseudopoten- 
tials (FA Ush and a kinetic energy cutoff of 30 Ry is used 
for the wavef unction expansion. Considering the fact that the 
force acting on H2 is much smaller than those in usual bonded 
systems, we set a very high total-energy convergence criteria, 
1.0 X 10~^ Ry, for the self consistent vdW-DF calculation to 
obtain accurate forces on atoms. The structures of the frag- 
ment models (with and without H2) are optimized by mini- 
mizing the forces on each atom to be smaller than 1.0 x 10~^ 
Ry/Bohr. The structures of the cry stall (with and without H2) 
are optimized by performing full variable-cell relaxations un- 
til the forces on atoms are smaller than 1.0 x 10~^ Ry/Bohr 
and the change in total energy of the primitive cell between 
two neighboring variable-cell steps is smaller than 1.0 x 10~^ 
Ry. A 4x4x4 Monkhorst-Pack k-point mesh is used to sam- 
ple the Brillouin zone of the crystall. 

III. RESULTS AND DISCUSSION 

A. Structure optimization with a single H2 adsorbed 

We first consider a single H2 molecule adsorbed on the frag- 
ment models and in the MOF-5 crystal structure. Since the 
physisorption is quite weak, the most stable adsorption con- 
figuration needs to be determined very carefully. For the frag- 
ment models we do the structure optimization in two steps. 
First, we perform a series of single point total-energy calcula- 
tions for different adsorption distances within a roughly cho- 
sen range (see Fig. O. Second, based on the lowest-energy 
adsorption distance obtained we then do a full structure re- 
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FIG. 3: (a) The total energy as a function of the adsorption distance 
between the center of the unrelaxed H2 and the Ce plane of the un- 
relaxed BDC model. The H2 and BDC model are fully optimized 
before the calculation, (b) Similar plotting for the corner model with 
a H2 adsorbed on the cup site. The adsorption distance is between 
the center of the H2 and the C3 plane. 



laxation with a very high convergence criteria, as mentioned 
previously, to ensure an accurate final adsorption configura- 
tion. The final optimized adsorption distances are also used 
as the initial input for the calculation of the crystal structure. 
For example, for the BDC model the result from the first step 
(see Fig.Oa)) gives a lowest-energy distance of 3.40A in line 
with the typical vdW interaction distance 1 16]. The final fully 
optimized distance from the second step is 3. 48 A. For the cor- 
ner model, the total-energy curve in Fig. Ob) shows a lowest- 
energy distance of 2.28A and the fully relaxed one is 2.0lA. 
The corresponding results for the crystal structure are 3. 49 A 
and 2. 06 A for the two adsorption sites, respectively. The full 
results are listed in Table HI 



B. Binding energy of a single H2 and the accuracy of the 
vdW-DF method 

We first examine the accuracy of the vdW-DF method by 
performing calculations for the fragment models and com- 
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corner+H2 BDC+H2 edge+H2 p(c)+H2 p(l)+H2 



distance 


2.01 


3.48 


3.46 


2.06 


3.49 


vdW-DF 


8.68 


4.55 


5.31 


9.39 


5.50 


CCSD(T) 




4.32^ 








MP2 


7.50" 


5.41^ 


5.10" 






PBE+D 


5.30" 






9.60" 





TABLE I: The calculated adsorption distance (distance, in A) 
and binding energy (vdW-DF, in kJ/mol) of a single H2 molecule 
adsorbed on the cup site or linker site of the fragment models 
(corner+H2, BDC+H2, and edge+H2) and crystal structure (p(c)+H2, 
p(l)+H2). Results of binding energy from some previous calculations 
(including CCSD(T), MP2, and PBE+D) are also Usted for a compar- 
ison. 

" data from Reference UM, 
^ data from Reference II29I] . 



pare the results of H2 binding energy to the existing results 
from the correlated quantum chemistry methods. The calcu- 
lated results are listed in Table. U together with reliable re- 
sults from CCSD(T), MP2, and DFT+D using PEE functional 
(PBE+D). Sagara et al has reported H2 binding energy for 
the BDC model using the CCSD(T) and MP2 methods adopt- 
ing the large QZVPP basis set |18] and got results of 4.32 
and 5.41 kJ/mol, respectively. Our vdW-DF result is 4.55 
kJ/mol which is very close to the CCSD(T) one. For the cor- 
ner model, a MP2 calculation Lllil using a basis set of aug- 
cc-pVQZ plus diffuse functions gave a result of 7.50 kJ/mol 
while a PBE+D calculation |TT] adopting a plane- wave basis 
set with 400 eV kinetic-energy cutoff showed 5.30 kJ/mol. 
Our result of 8.68 kJ/mol is closer to the MP2 result than 
the PBE+D one does. For the edge model, the present re- 
sult of 5.31 kJ/mol is also quite close to the MP2 result 5.1 
kJ/mol obtained in Ref. illll . From the above comparison, 
we can conclude that the vdW-DF method can describe well 
the physisorption of H2 molecule in MOF-5 and its accuracy 
is comparable to that of the MP2 and CCSD(T) methods and 
seems better than the DFT+D. An outstanding advantage of 
the DFT-based methods is that, unlike the correlated quantum 
chemistry methods, they are not so strongly dependent on the 
basis set used. For the latter, although the theories themselves 
are more accurate but their results are very scattered in the 
literature, varying largely with different basis sets used, and 
only become reliable when the very big basis sets (as men- 
tioned above) are adopted. Consequently, even for the small 
fragment modes they are computationaly very demanding. 



C. Fragment models vs. crystal structure 

As one can see in Table II the fragment modes and the 
crystal structure give quite different results of binding energy 
for a single H2 adsorbed. The binding energy with the cor- 
ner model (corner+H2) is 0.71 kJ/mol lower than that with 
the crystal structure (p(c)+H2) and the binding energy with 
the BDC model (BDC+H2) is 0.95 kJ/mol smaller than the 
crystal one (p(l)+H2). This difference is due to the effect 
from the extended environment in the crystal structure which 



may provide an accumulated electrostatic and vdW forces be- 
cause of their long-range nature. Additionally, there will be 
some charge transfer between a fragment and its environment, 
which will also change the strength of the physisorption. Our 
result shows that this effect is particular significant when the 
H2 molecule is adsorbed on the organic linker part: Without 
connecting to the metal oxide corner part the BDC model un- 
derestimates the binding energy by about 20%. After this con- 
nection is made, as in the case of the edge model (edge+H2), 
then the result becomes closer to (but still smaller than) that 
of p(l)+H2. However, this connection will increase signifi- 
cantly the size of the model and make the higher-order quan- 
tum chemistry CCSD(T) calculation a big challenge. 

Our calculation demonstrates that, compared to the true 
crystall structure calculation, the small fragment modes can 
still capture the main feature of the physisorption of H2 on 
the different adsorption sites: Both show that the binding on 
the cup site of the 0Zn4 corner is much stronger than that on 
the linker site. This has also been comfirmed experimentally 
by inelastic neutron scattering spectroscopy ||3Q|i and neutron 
powder diffraction measurement f37]. The physics underly- 
ing is that the metal ion induced charge polarization around 
the corner causes stronger electrostatic force on the H2. How- 
ever, our calculation shows that the small fragment models ig- 
noring the influence of the adjacent components are not quan- 
titatively sufficient to represent the different parts of MOFs 
for investigating the hydrogen storage with MOFs. To get ac- 
curate H2 binding energies for the different adsorption sites, 
calculations with the true periodic boundary conditions for the 
crystal structure are needed. Besides the underestimated bind- 
ing energy for a single H2 adsorbed, the fragment models also 
give a different trend for multiple H2 adsorbed, as will be dis- 
cussed later. 



D. Effect of multiple H2 adsorbed 

In the practice of hydrogen storage many H2 molecules will 
be adsorbed in the primitive cell of a MOF crystall. Here, we 
investigate the effects of hydrogen density on the binding en- 
ergy by considering multiple H2 molecules adsorbed in the 
MOF-5 crystall, on the corner site (p(c)) and the linker site 
(p(l)), as shown in Fig.|4l The lattice constant and the average 
binding energy per H2 for difference H2 densities are listed in 
Table [III One can see that when the number of H2 adsorbed 
in the primitive cell increases, the lattice constant actually de- 
creases a little bit (by 0.02 A) and the binding energy also de- 
creases by 0.34 kJ/mol and 0.36 kJ/mol for the p(c) and p(l), 
respectively. This small decrease of the lattice constant im- 
plies that there exists an attractive force among the adsorbed 
H2 molecules, which increases with the increaseing H2 den- 
sity. This behavior is different from the case of chemisorp- 
tion where an increase in adsorbate density will usually lead 
to an expansion of the volume. For the highest H2 density 
considered in the present work, the smallest H2...H2 distance 
is around 4.1 A which is quite close to that in the bulk H2 
solid (3.8 A) 132]. For even higher H2 density comparable 
to that in the solid state, a previous calculation by Yildirim et 
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FIG. 4: The MOF-5 primitive cell with (a) one H2 adsorbed on the cup site, (b) 8 H2 adsorbed on the cup sites, (c) one H2 adsorbed on the 
linker site, and (d) 12 H2 adsorbed on the linker sites. 





corner+H2 




corner+4H2 


BDC+H2 




BDC+2H2 


distance 


2.05 




1.93 


3.48 




3.46 


vdW-DF 


8.68 




8.78 


4.55 




4.56 




P(C)+H2 


F-p(c)+H2 


p(c)+8H2 


P(l)+H2 


F-p(l)+H2 


p(l)+12H2 


distance 


2.06 


2.05 


2.07 


3.49 


3.45 


3.48 


lattice const. 


18.81 


18.73 


18.79 


18.81 


18.72 


18.79 


vdW-DF 


9.39 


9.82 


9.05 


5.50 


10.39 


5.14 



TABLE II: Calculated average binding energy (vdW-DF, in kJ/mol) per H2 for one and multiple H2 adsorbed on the fragment models (the first 
part of the table) and in the cry stall structure (the second part of the table). The corresponding adsorption distance and the lattice constant are 
also listed. F-p(c)+H2 (F-p(l)+H2) denotes one H2 adsorbed on the cup (linker) site in the cry stall structure with the fluorine substitution, as 
described in the text. 



al showed that a unique three-dimensional interlinked H2 
nanocage will be formed, which will enhance further the in- 
teraction among the H2 molecules. Since the attraction among 
the H2 molecules is in the opposite direction to their binding 
to the MOF-5 crystall, it reduces the binding energy conse- 
quently. 

As a comparison to the crystall case, we also check the 
multiple H2 effect using the fragment models. As shown in 
Table Hn for the corner model, when the four cup sites are oc- 
cupied by H2 , the average binding energy actually increases 
by 0.1 kJ/mol compared to the case of only one H2 adsorbed. 
For the BDC model, this change is much smaller but is still 



increased (by 0.01 kJ/mol) with H2 adaorbed on both sides of 
the benzene ring compared to only on one side. This quali- 
tative difference between the crystall and the fragment mod- 
els for multiple H2 adsorption can be understood by consid- 
ering the fact that in the fragment models the H2 molecules 
are adsorbed only on the surface of the cluster, and therefore, 
the interlink between H2 molecules adsorbed in a cage of the 
crystall structure is absent. In this case, when multiple H2 
molecules are adsorbed the synergistic effect of electrostatic 
and vdW forces leads to the small increase of the binding en- 
ergy. Thus, our calculation shows that, compared to the true 
crystall stucture, the fragment models not only underestimates 



FIG. 5: Primitive cell with the H atoms of the organic linkers replaced by F atoms: (a) with one H2 adsorbed on the cup site, (b) with one H2 
absorbed on the linker site. 



the binding energy when a single H2 molecule is adsorbed but 
also give an opposite trend of the binding energy when multi- 
ple H2 molecules are adsorbed. 

Another important effect of multiple H2 adsorption is an 
perpendicular-to-parallel rotation of the H2 molecules ad- 
sorbed on the cup site: When only one H2 is adsorbed, it 
prefers to be perpendicular to the C3 plane (see Fig. lUa)). 
However, when multiple H2 molecules are adsorbed on 
the cup site, their adsorption configuration will undergo an 
perpendicular-to-parallel transition and becomes parallel to 
the C3 plane (see Fig.Ub) ). This transition reflects the inter- 
action between the H2 molecules adsorbed on the same corner 
unit. 



E. Substitution of F atoms for the H atoms of the organic 
linker 

As is already shown, the H2 binding energy with the or- 
ganic linker is significantly smaller than with the metal ox- 
ide corner, which largely limits the H2 loading for hydrogen 
storage application. The reason is that the electrostatic force 
coming from the 0Zn4 cluster is much stronger than the vdW 
force coming from the organic linker. This implies a way for 
improving the average binding energy of H2 molecules in the 
MOF-5 cry stall: It may be achieved by creating extra polar- 
ization in the organic linker part and therefore enhancing the 
binding energy by the induced electrostatic force. 

Here we consider a possibility that the 4 hydrogen atoms 
of an organic linker are replaced with 4 fluorine atoms (see 
Fig. [5]) which have a much stronger electronegativity leading 
to a larger charge transfer and, therefore, will give rise to ex- 
tra polarization. Fluorine atom has also a small mass for a 
good gravimetric hydrogen loading. The results of binding 
energy after the fluorine substitution are also given in Table. [Ill 
(F-p(c)+H2 and F-p(l)+H2). One can see that for one H2 ad- 
sorbed on the cup site the binding energy increases slightly 
from 9.39 kJ/mol to 9.82 kJ/mol. However, for one H2 ad- 
sorbed on the linker site the binding energy is enhanced sig- 



nificantly from 5.50 kJ/mol to 10.39 kJ/mol which is now even 
larger than the binding energy with the corner, showing that 
the proposed extra electrostatic interaction is really working. 
This result casts a light on possible ways for improving the 
performance of MOFs for hydrogen storage. 



IV. SUMMARY 

By using the vdW-DF method implemented with plane- 
wave basis functions we have investigated systematically the 
physisorption of hydrogen molecules in MOF-5 crystall as 
used for hydrogen storage. We first examed the accuracy of 
this methodology by comparing the results for three framgent 
models to those from the correlated quantum chemistry meth- 
ods, CCSD(T) and MP2, as well as from the DFT+D method. 
Good accuracy was found, which is comparable to that of the 
CCSD(T) and MP2 with very large basis sets. By compar- 
ing the results for the fragment modes with those for the true 
crystall structure we have shown that the fragment models un- 
derestimate largely the binding energy and therefore cannot 
represent well the hydrogen- storage property of the crystall 
which, however, cannot be dealt with by the quantum chem- 
istry methods. For many hydrogen molecules adsorbed, the 
fragment models even give an opposite trend of the binding 
energy. 

Our calculations for the true crystall structure show that 
the binding energy with the metal oxide corner is signifi- 
cantly larger than with the organic linker due to the metal ion- 
induced electrostatic force. With many hydrogen molecules 
adsorbed the crystall volume and the average binding energy 
decrease a little bit because of the small attraction among the 
hydrogen molecules, which leads to the regular interlinked H2 
nanocage structure for higher H2 density found previously. 
For many hydrogen molecules adsorbed on the corner they 
undergo a perpendicular-to-parallel transition and prefer to be 
parallel to the C3 plane. 

The present work shows that the much weaker binding with 
the organic linker can be improved significantly by introduc- 
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ing extra polarization and the resulting electrostatic interac- 
tion: The substitution of F atoms for the H atoms can increase 
the binding energy from 5.50 kJ/mol to 10.39 kJ/mol which is 
now very close to the binding energy with the corner. This 
casts a light on possible ways for improving the hydrogen 
loading of MOFs for hydrogen storage. 
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